Bar-mode instability of rapidly spinning black hole in higher dimensions: 
Numerical simulation in general relativity 
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Numerical-relativity simulation is performed for rapidly spinning black holes (BHs) in a higher- 
dimensional spacetime of special symmetries for the dimensionality 6 < d < 8. We find that higher- 
dimensional BHs, spinning rapidly enough, are dynamically unstable against nonaxisymmetric bar- 
mode deformation and spontaneously emit gravitational waves, irrespective of d as in the case d = 5 
The critical values of a nondimensional spin parameter for the onset of the instability are 
q ~ a/ fi^^^'^''^'' « 0.74 for d = 6, ~ 0.73 for d = 7, and ~ 0.77 for d = 8 where /i and a are mass 
and spin parameters. Black holes with a spin smaller than these critical values (gcrit) appear to 
be dynamically stable for any perturbation. Longterm simulations for the unstable BHs are also 
performed for d = 6 and 7. We find that they spin down as a result of gravitational-wave emission 
and subsequently settle to a stable stationary BH of a spin smaller than qait- For more rapidly 
spinning unstable BHs, the timescale, for which the new state is reached, is shorter and fraction of 
the spin-down is larger. Our findings imply that a highly rapidly spinning BH with q > Qcrit cannot 
be a stationary product in the particle accelerators, even if it would be formed as a consequence of 
a TeV-gravity hypothesis. Its implications for the phenomenology of a mini BH are discussed. 

PACS numbers: 04.25. D-, 04.30.-w, 04.40.Dg 
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I. INTRODUCTION 

Clarifying formation and evolution processes of mini 
black hole (BH) in higher-dimensional spacetimes has be- 
come an important issue since a possibility of BH forma- 
tion in huge particle accelerators was pointed out. If our 
space is a 3-brane in a higher-dimensional flat spacetime 
of spacetime dimensionality d > 6 or in an anti-de 
Sitter (AdS) spacetime of d > 5 the Planck energy 
could be of O(TeV) that may be accessible with parti- 
cle accelerators in operation, the CERN Large Hadron 
Collider (LHC). In the presence of the extra dimensions, 
BHs of very small mass energy > TeV may be produced 
during the particle collision in the accelerators because 
the true Planck energy may be as low as TeV scale. 

A hypothetical phenomenology of a BH produced in 
the huge particle accelerator was first discussed in [3, [1| 
(see [11 for reviews). According to this standard scenario, 
a mini BH evolves in the following manner: During the 
high-energy particle collision of a sufficiently small im- 
pact parameter and of energy sufficiently higher than 
the Planck energy, two particles will merge to form a 
deformed BH, and then, it settles to a quasistationary 
state after emission of gravitational waves. The typi- 
cally assumed timescale for gravitational-wave emission 
is ~ 10r_|_ /c (about 10 times of the dynamical timescale) 
where r_|_ and c are the horizon radius and speed of 
light, respectively. The quasistationary BH will be subse- 
quently evaporated by the Hawking radiation Q , because 
of quantum-field-effects in a curved spacetime. Much ef- 
fort has been devoted to calculating the greybody factor 
in the Hawking radiation for improving the prediction 
of si gna ls in the particle-collision experiments |8l-[l3| (see 
also for related issues). 

By contrast, the analyses for BH formation after the 



particle collision and subsequent evolution by gravita- 
tional radiation reaction in higher-dimensional spacetime 
have not been done yet [but see [Tsl - flU for studies in 
the four-dimensional (4D) case]. These phases are ex- 
pected to be described well in the context of general 
relativity ^13], but due to its highly nonlinear nature, 
any approximation breaks down. Obviously, numerical- 
relativity simulation is the unique approach for studying 
this phase. 

One of the important issues has been to clarify what 
type of BH is formed and whether it is stable or not. In 
the 4D case, any stationary BH formed in vacuum has to 
be a Kerr BH (neglecting the electric ch arge of the BH) 
because of the uniqueness theorem (e.g., [20| for review), 
and the Kerr BH has been proven to be stable [2l| - [23| 
(but see [l^l for remaining issues for a perfect proof). 
These facts strongly constrain the possible scenario for 
mini BH formation and its subsequent evolution. By 
contrast, there is no uniqueness theorem and no proof 
for the stability of higher-dimensional BHs (but see [1^ 
for uniqueness of 5D BHs of the spherical horizon topol- 
ogy). As a result, the standard scenario described above 
(i.e., formation, evolution by subsequent gravitational- 
wave emission, and evaporation by the Hawking radia- 
tion) is quite uncertain. 

A mini BH, if it is formed as a result particle colli- 
sion in a higher-dimensional spacetime, will have only 
one spin parameter associated with the orbital motion. 
This can restrict the possibility for the type of the formed 
BH. However, even in this case, one cannot restrict the 
BH type in higher dimensions. For example, in the 5D 
case, there are many types of BHs: e.g., usual Kerr-type 
BH (Myers-Perry BH [26|) for which surface of the event 
horizon is S^, fat and thin black rings for which surface 
has a ring-like shape [l^ . Recently, black di-ring [2^ and 
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black-saturn [29| solutions were derived as well (see also 
(3C| for other exact solutions of 5D black objects but with 
two spin parameters) . Several authors 0, [sil, [13] discuss 
a possibility of the black-ring formation in particle colli- 
sions. Although it is still an open question, the analysis 
of apparent horizons in 31] indicates that black-ring for- 
mation is not very likely in two point-particle system. 

Higher-dimensional BHs with a high spin parameter 
are known to be unstable against axisymmetric pertur- 
bations. Eniparan and Myers [s^ speculated that rapidly 
spinning BHs with the spacetime dimensionalit y d > 6 
are subject to the Gregory-Lafiamme instability [34| , be- 
cause they have a high degree of ellipticity (i.e., the black 
membrane limit). Very recently, Dias et al. and Murata 
et al. indeed showed, by a linear perturbation analy- 
sis, that rapidly spinning BHs for 6 < d < 9 are unsta- 
ble against axisymmetric multiple-ring-like deformation 
[H,!!!] (see also related papers [STj). 

On the other hand, little is known for the stability 
of spinning BHs against nonaxisymmetric perturbations. 
Emparan and Myers [SS*] also speculated, based on a ther- 
modynamic argument, that rapidly spinning BHs with a 
sufficiently high spin may be unstable against nonaxisym- 
metric perturbation for d > 5, because the horizon area 
(the so-called entropy) of a rapidly spinning BH is often 
smaller than that of two boosted Schwarzschild BHs with 
the same total energy and angular momentum. How- 
ever, the correspondence between the thermodynamical 
and dynamical instabilities has not been well established. 
Authors in (Hj guess, based on a hydrodynamic/gravity 
correspondence argument, that a nonaxisymmetric insta- 
bility may occur for spinning BHs with spin smaller than 
the critical value for the onset of axisymmetric insta- 
bilities. However, a rigorous and quantitative analysis 
is absent in these studies. To strictly clarify what the 
criterion for the onset of nonaxisymmetric instabilities 
is and how the instabilities occur and proceed, we have 
to solve Einstein's equation. Recently, we performed a 
numerical-relativity simulation for 5D spinning BHs for 
the first time, and found that BHs spinning rapidly are 
dynamically unstable against nonaxisymmetric bar-mode 
deformation if the spin parameter satisfies the condition 
q :— a/ 11^^^ > 0.87 (see Sec. H for the definitions of q, a, 
and n) Jj]: An interesting fact is that the critical value 
we found is close to the value predicted by Emparan and 
Myers |33|, ss 0.85, suggesting that their argument rely- 
ing on BH thermodynamics may be reliable (although we 
show in this paper that this is not the case for d > 6). 

In this paper, we study the stability of rapidly spin- 
ning BHs with dimensionality d > 6 using a new 
numerical- relativity code SACRA-ND, which is extended 
from SACRA ^ and SACRA5D 1]. In this code, Ein- 
stein's equation for higher-dimensional spacetimes of any 
dimensionality for d > 5 is fully solved without impos- 
ing axial symmetry. Thus, in this work, we do not have 
to assume that the amplitude of the perturbation from 
the background axisymmetric BH solution is small nor 
restrict attention to an axisymmetric perturbation. The 



merits of this approach are that (i) a wide variety of insta- 
bilities can be investigated with no ambiguity and with 
no approximation (except for finite-difference approxima- 
tion); (ii) the final fate after the onset of the instabilities 
can be determined because the amplitude of the pertur- 
bation from the background BH solution does not have 
to be small; and (iii) the criterion for the onset of the 
instabilities is quantitatively determined. In this paper, 
we focus on the six-, seven-, and eight-dimensional (6D, 
7D, and 8D) BHs of single spin parameter because such a 
class of BH is a possible outcome in the particle acceler- 
ators. We shall show, for the first time, that the rapidly 
spinning BHs are dynamically unstable against bar-mode 
deformation irrespective of the dimensionality. We, fur- 
thermore, evolve the 6D and 7D unstable BHs for a long 
time until they settle to a new stable state. This enables 
to quantitatively clarify a possible (classical) evolution 
process of an unstable mini BH for the first time. 

The paper is organized as follows. In Sec. II, the formu- 
lation for solving Einstein's evolution equation in higher 
dimensions is described. In Sec. Ill, a method for evolv- 
ing Einstein's equation in the spacetime of a special sym- 
metry is summarized. In Sec. IV, we present the numer- 
ical results focusing on the criterion for the onset of the 
bar-mode dynamical instability of rapidly spinning BHs 
and the fate after the onset of the instability. We will 
show that the rapidly spinning BHs, which have a spin 
larger than the critical value, are unstable against spon- 
taneous emission of quadrupole gravitational waves. We 
will also show that all the unstable BHs, considered in 
this paper, evolve as a result of gravitational radiation 
reaction, settling to a stable BH of spin smaller than the 
critical value. In Sec. V, the conditions for the sponta- 
neous gravitational-wave emission and associated spin- 
down process are clarified. Section VI is devoted to a 
summary and discussion. In particular, implications of 
our results for the phenomenology of a mini BH are dis- 
cussed in detail. 

In Sees. I - VI B, we adopt the units in which c = 1. 
Only in Sec. VI C, the natural units c = ^ = 1, where H 
is the Planck constant, are adopted, d denotes the space- 
time dimension with n := d — A being the number of the 
extra dimension, and Gd denotes the d-dimensional gravi- 
tational constant. The Cartesian coordinates {x, y, z, Wq) 
are used for the space coordinates: = {x,y, z) denote 
the usual three dimensional coordinates, and Wq {q = 1 
- n) the coordinates of extra dimension, t denotes the 
time. Indices i, j, k, and I denote the general spatial 
coordinates. 



II. FORMULATION 

We consider a vacuum higher-dimensional spacetime of 
SO(n-|- 1) symmetry for which the line element is written 
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(1) 



where a is the lapse function, /3'^ the shift vector, the 
space metric, and dfl^ the Une element of n-dimensional 
unit sphere. 7„„ is a conformal factor for the extra- 
dimensional metric components. x° denotes (x, y, p) 
where p is a radial coordinate 



9=1 



(2) 



The geometric quantities in this symmetric spacetime de- 
pend only on t, x, y, and p. (Note that by SO(n -I- 1) 
symmetry, we imply that the subspace of {z,wi, - ^Wn) 
coordinate directions is isotropic.) 

The line element of a Myers-Perry BH of single spin 
(i.e., the Kerr- type BH) in the Boyer-Lindquist-type co- 
ordinates is USi 



ds^ = -dt^ 



,, ^-(dt — a sin Odw) + -rdf 



+ (f2 + a^) sin^ edip^ + f2 cos^ BdUl, 



(3) 



where p and a are mass and spin parameters, respec- 
tively, E := + cos^ 6*, and A :— P -\- — pf^^^. 
Note that the mass and angular momentum of this BH 
are 



M = 



J 



{d-2)nd-2P 



d-2 



Ma, 



(4) 
(5) 



where is the area of (d— 2)-dimensional unit sphere, 
i.e., = 27r('^-i)/Vr[(d- l)/2]. Thus, the line ele- 

ment denoted by Eq. ([T]) includes the spinning BH solu- 
tion (|3]), and nonstationary, nonaxisymmetric deformed 
states can be described as well. 

In numerical simulation, we adopt the Cartesian coor- 
dinates {x, y, z, Wq) instead of the curvilinear coordinates; 
e.g., for d = 7, the relations between (z, wi, u'2, W3) and 
radial and angular coordinates (p,ip,tpi,ip2) are 



z — p cos tp, 

wi — psinip cosifii, 

W2 = p sin sin ipi cos ip2 , 

W3 ^ p sin ip sin tpi sin ip2 ■ 



(6) 
(7) 
(8) 
(9) 



Then, Einstein's evolution equation is solved in the 
Cartesian coordinates using the so-called cartoon method 
[40, 43]. Namely, we solve the equations in the (x,y,z) 
hyperplane; the hyperplane of ip ~ — ^- The method 
used in the present work will be described in the next 
section. 

To solve Einstein's evolution equation in the Carte- 
sian coordinates, we adopt a multidimensional version of 



the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) for- 
malism (4G . 1411] (see also [4§| for a different formalism) . 
We rewrite the line element in the form, 

ds^ = ~{a^ - PkP^)df + 2Pkdx''dt + x^^lijdx'dx^ {IQ) 

where x = [det(7ij)]^^/('*^^) is a conformal factor and 
the conformal spatial metric 7ij satisfies the condition 
det(7.y ) = 1. In addition, we define the following quan- 
tities from the extrinsic curvature Kij, 



1 



as well as an auxiliary variable 

P := -djfi. 



(11) 
(12) 

(13) 



Then, the variables (x, 7jj , K, ^4^ , F*) are evolved solving 
the following equations [401 : 



(dt 
(dt 



d-1 



l3^dk)K = -D.D'a + a A'^ A 



K' 



d-l 



(14) 

(15) 
(16) 



( 



KA,, 



2A,kA' 



+ A,kd,p'' + Ak,d,p 



1 



(17) 



d~ 1 

2a 



J- jk^ 



~2A'^d^a-Pdj/3' 
-^^f''dkd,l3^+j^''d,dkf3' 
d-2,.._. {d~l)d,x. 



d-l 



(18) 



Here, Di and Rij denote the covariant derivative and the 
Ricci tensor with respect to 7^^, and TF denotes taking 
the trace-free part. 

For evaluating RjP^ , we first decompose Rij as 



Rij ^ij 



R^ 



(19) 



where Rij is the Ricci tensor with respect to 7ij and we 
write it in the following form 



R^3 = -^l'"l^]M + \ {lk^^JT'' + ^jd^V 



-pi -pk 

'^ik'^jl- 



(20) 
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Here, f^j, is the Christoffel symbol with respect to jij. 
Rf, is sum of the terms associated with y, 



(d-3) (Ax) Ax 



2X 

Then, we obtain 



2x 



4x 



^(Ax)^'x 



(21) 



3 

d- 3 
4x^ 



1 



d- 1 

b^b,x 



d- 



(Ax) Ax - ^7., (Ax)^'x 



.(22) 



For the 4D case, the so-called puncture gauge con- 
ditions are known to be robust for evolving BH space- 
time with the BSSN formalism This is also the case 
for the higher-dimensional spacetime [13] ■ However, the 
freely chosen coefficients in this gauge have to be care- 
fully determined for the stable and longterm evolution. 
Specifically, we choose the equations in the form 



dta = —l.baK, 



or 



(23) 
(24) 
(25) 



(26) 



where i?* is an auxiliary function, and rjB is a nondi- 
mensional constant for which we give different values for 
different number of d because a small value of t^b is not 
allowed for a high number of d or for a high spin for 
achieving the longterm stable numerical evolution. For 
d = 5, we gave ry^ = 1 in the previous work In the 
present work, we employ rys = 2 - 5 for d = 6, 3 - 5 for 
d = 7, and 8 for d — 8. For unstable BHs, a large value of 
r]B is favored in performing a longterm simulation until 
the growth of the instability saturates and subsequently 
the deformation damps. We tried to use both Eqs. (|25)) 
and , and found that both of them work well as far 
as BH spin is not extremely large. For a very high spin 
(g > 1), however, Eq. (|26p works better than Eq. (1^51) for 
a longterm stable simulation. 



III. MODIFIED CARTOON METHOD 

We solve Einstein's evolution equation in the space- 
time of SO(n -|- 1) symmetry with n = d — 4 using 
the Cartesian coordinates {x,y, z,wi, ■ ■ ■ ,w„) with the 
BSSN formalism. The SO(n -t- 1) symmetry is imposed 
for the {z,wi,--- ,w„) subspace (i.e., the subspace of 



(z,wi, ...jWn) coordinate directions is assumed to be 
isotropic). The method is qualitatively the same as those 
described for 5D spacetime [13]. Namely, we solve the 
equations in the {x,y,z) hyperplane (i.e., wi — W2 — 
■■■ = Wn = hyperplane), and derivatives with re- 
spect to these coordinates are evaluated by a straightfor- 
ward finite differencing. On the other hand, the deriva- 
tives with respect to the extra-dimensional coordinates, 
Wq {q — l---n), are evaluated using symmetry rela- 
tions. In the previous paper [40|, we adopted the original 
prescription in the cartoon method often employed for 
axisymmetric spacetimes [43|: We prepare 4 additional 
grid points for the fifth coordinate, w, i.e., 2 ± Aw and 
zLAw. Then, the values for all the geometric quantities 
at these grid points are determined using the symmetry 
relation associated with the Killing vector d/dip where 
-0 ~ tan^^(w/z), and then, the derivatives with respect 
to w arc evaluated using a fourth-order finite differencing. 

This original method is quite simple for code imple- 
mentation. However, it is memory-consuming in higher- 
dimensional simulations, because we have to increase the 
grid number by a factor of five whenever the number of 
extra dimension is increased. To save the memory used in 
computation, we employ a different method in which the 
derivatives with respect to Wq are replaced to those with 
respect to z without preparing additional grid points for 
the extra-dimensional direction. In the following, we de- 
scribe the new method denoting the scalar, vector, and 
tensor by Q, Q\ and Qij, respectively, and decomposing 
the subscripts i,j into a, b = x, y, z and Wq {q = I - n). 

In the (a;, y, z) hyperplane, the following relations hold 
because of the assumed isotropy: 



and 



(27) 



(28) 



The derivatives with respect to a;" of Eq. ([28)1 are trivially 
zero. For other derivatives, the following relations hold: 
For the scalar quantities. 



{Q,z/z)&qr 



(29) 



for the vector quantities, 

Q ,Wq ~ Q ,bWq ~ Q ''.WqWr 
Q ,WqWr ~ {Q -zI ^ ^1^^ 

Q%v-. = (^?V^),.V. (30) 

and for the tensor quantities, 

Qab,Wq — QwpWq — 0, Q Awg^w^ — Az / ^')^qr^ 

Q ZWq^Wr [{QzZ Qwti) ) / ^] ^(jr ; 

QAB.WqWq — QaB.z/ QAz.WqWq — {QAz/z),Z-, 
Q ZZ^WqWq QzZ^z/ ^ ~^ (^/^ }(^QzZ ~^ Q Ww) 1 

Qww .z I ^ ^^ (2/ Z '){C^zz ~^ Qiuw): 

Qww,z/z (q^r). (31) 



fWqWq.WqWq 
^WqWa-Wj-Wr 
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Here, A, B = x or y (not z) and Sgr is the Kronecker's 
delta. We did not evaluate Qij.aw, and Qij^wqw,- (q 7^ 
r) even if they are not vanishing, because they do not 
appear in the BSSN equations with SO(n + 1) symmetry: 
Note that the second derivatives appear only in the term 
-{l/2)^'''-dkdi% of Eq. ([201) and 7^"' satisfies the relation 



of Eq. 



. We also do not have to evaluate Q 



aWq ,ZJ 



and 



QwgWr,ij (q 7^ because they appear only in the aWq 
and WqWr components of Eq. (ITTl) with Eq. (pO)) which 
do not have to be evolved, as mentioned later. 

Some of the above prescriptions can be used only for 
z 7^ in computer because of the presence of the terms 
associated with 1/z, although these terms are actually 
regular. Thus for 2 = 0, the following relations, which 
are found from SO(n + 1) symmetry, are employed: 



0, 



Q %q ~ Q ,Z-) Q ,WqWq 

Q AB ,WqWq Q AB ,ZZ ■) Q AZjWqWq O7 

Q AWq,Wq QAZjZj Q ZWq.Wq ^ 0; 

Q ZZjWqWq QwW^ZZj QwqWq.WqWq QzZ,ZZl 

QwqWq,Wr-Wr QwqWq,ZZ- 

(32) 



The final remark is on the treatment of the advection 
terms such as fi'^dtjij- Because (3"^" — for the Wq — 
hyperplane in the present case, the advection terms asso- 
ciated with is always vanishing in the computational 
domain chosen in our method. 

With these prescriptions, all the derivatives associated 
with the extra-dimensional coordinates can be replaced 
to the finite-differencing terms with respect to or with 
no finite differencing. It is worthy to note that total 
amount of computational operation is only slightly larger 
than that for the 3-1-1 case. 

Because of the relations (P7|) and I^E^ that follow from 
SO(n-|- 1) symmetry, the implementation for the higher- 
dimensional contribution can be even simplified: We have 
only to evolve the scalar equations, the "a" components 
of the vector equations, and the "ab" components and 
(one of) the WqWq components of the tensor equations 
(i.e., the equations for and Aij) of the BSSN formal- 
ism. Therefore, we need to increase only one component 
for 7ij and Aij irrespective of the dimensionality. In the 
evolution equation, we often have terms such as /3* ^ or 

I'^^lak,!, which are evaluated by 



(33) 



where the prescriptions shown in Eqs. - l\'3'2\i are used 
for evaluating the second terms in the right-hand side. 
These facts imply that once a 5D code is implemented, it 
is quite straightforward to extend it to a code for d > 6, 
even when we do not employ a curvilinear coordinate 
system [42j . 



IV. NUMERICAL SIMULATION 

A. Setting and Methodology 

We prepare 6D - 8D Myers-Perry BHs of single spin 
parameter (see Eq. ([3]) for the line element in the Boyer- 
Lindquist coordinates) as the initial condition. Because 
the Boyer-Lindquist coordinates are not suitable for the 
Cauchy evolution of BH spacetimes, we transform the 
radial coordinate introducing a quasiradial coordinate r; 
we rewrite the line element in the form (see |45j] for the 
4D case) 



ds' = -dt' 



„2 I „2 i/i2 



(dt - asin^ 9d(pf 



+ <^>{dr^ +r^d9^ 
+ (f2 + a^) sin^ 6dip^ + cos^ 9dnl, (35) 

where $ is a conformal factor for the (r, 6) plane. Namely, 
the two-dimensional metric for (r, 0) is written in a con- 
formally flat form, and the following relations are satis- 
fled: 



$i/2dr ==±(E/A)i/2df, 
$1/2^ ^ 



Then, r is defined by 



r^^exp 



± 



dR 



-d 



(36) 
(37) 



(38) 



where the plus and minus signs are adopted for the re- 
gions r > rfi and r < rh, respectively. r+ is the horizon 
radius in the Boyer-Lindquist coordinates and a positive 
root of -\- a? = jir^'^ . is the horizon radius in 
the quasiradial coordinate, which is determined by the 
condition r = f for f — > 00. A fourth-order numerical 
integration is performed for the integral of Eq. ([55]) us- 
ing the Bode's rule ^] because it cannot be analytically 
integrated in general. 

With this transformation, the t =const hypersurface 
becomes spacelike everywhere for < r < 00, and fur- 
thermore, the singularity, which is originally located at 
f = 0, is excluded. More specifically, the spacelike hy- 
persurface has an inversion symmetry with respect to 
r = Th hypersphere (i.e., this hypersphere is the worm- 
hole throat), and thus the "point", r = 0, represents 
the spacelike infinity of (say) another world beyond the 
horizon. Although this point is not a physical singu- 
larity, it becomes a coordinate singularity, because $ 
is proportional to irrespective of the number of d. 
However, the puncture approach [specifically, appropri- 
ate choice of the conformal factor in numerical simula- 
tion, X — [det(7y )]^^/('*^^\ and choice of the puncture 
gauge ((23|) - ((25)) ] enables to stably evolve a BH space- 
time with no difficulty [13, El] . 

The condition for the onset of instabilities for spin- 
ning BHs depends on a nondimensional spin parameter. 
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Because the stability does not depend on the magni- 
tude of /X, all dimensional quantities can be scaled out 
of the problem appropriately normalizing them by using 
/i: Specifically, j^g^g ^j^g dimension of length and 

time in the c — 1 units (and is often referred to as Ts, 
e.g.,[ll[), and thus, we should define a nondimensional 
spin as 



q - 



l/{d-3) ■ 



(39) 



The stability of a Myers-Perry BH depends only on this 
quantity. Note that this is different from a* := q/r+ 
which is often used as another nondimensional spin pa- 
rameter 0, [ll[ . 

Numerical simulation is performed in the {x, y, z) co- 
ordinates, and thus, the initial condition is prepared by 
performing a coordinate transformation from (r, 0, ip) to 
{x,y^z). In addition, we add a small bar-mode pertur- 
bation to the conformal factor x a-s 



1 + A- 



2/{d-3) \ 2rl 



(40) 



where xo is the nonperturbed solution and A the initial 
perturbation amplitude, A <^ 1. We focus here on the 
bar mode because it is often the most relevant unstable 
mode for self-gravitating, dynamically unstable rotating 
systems such as a rotating star (47|- We chose the value 
of A for a wide range from 10^^ to 0.02 and found that for 
unstable BHs, the growth rate of the unstable mode does 
not depend on the initial magnitude of A. Taking into 
account this result, in a longterm evolution of unstable 
BHs for studying the nonlinear growth of the unstable 
mode, we choose a relatively large value oi A = 0.005 
or 0.02 to save computational costs. It should be also 
pointed out that the numerical error accumulates with 
time and a significant resolution is needed to suppress 
it. This implies that if a simulation was started with 
a very small value of A, quite expensive computational 
costs would be required to accurately follow the nonlinear 
growth of the perturbation. 

Initially, the lapse function is chosen as a = x- With 
this modification, a becomes positive except for r = 
where a = (near r = 0, a is proportional to with 
this choice). On the other hand, the shift vector is not 
modified; we choose the same shift as in the Myers-Perry 
BH initially. 

We solve Einstein's evolution equation in the BSSN 
formalism using a new code, SACRA-ND, in which an 
adaptive mesh refinement (AMR) algorithm is imple- 
mented in the same manner as in SACRA [39]. Be- 
cause we solve the higher-dimensional equations only 
with the {x, y, z) coordinates, the methods of interpola- 
tion, extrapolation, and evolution for the grid structure 
in SACRA-ND are totally the same as those in SACRA: 
All the spatial derivatives with respect to = {x, y, z) 
are evaluated using a centered fourth-order finite differ- 
encing except for the advection term such as P^dkX for 



which a fourth-order upwind scheme is adopted. The 
time evolution is carried out using the standard fourth- 
order Runge-Kutta method The reader may refer 
to (39| for details about the numerical methods. 

For the AMR scheme, we prepare 6 refinement levels 
in the present numerical simulation. As a test, we per- 
formed several simulations enlarging the computational 
domain with 7 refinement levels while the grid resolution 
in the finest level is unchanged. We found that the re- 
sults shown in Sec. IV are essentially independent of the 
location of the outer boundary. Computational domain 
for each refinement level has a half cubic shape which 
covers [—Li : Li] for x and y, and [0 : Li] for z; we as- 
sume the equatorial plane symmetry. Here, Li denotes 
the location of the refinement boundary for the l-th re- 
finement level with I = - Ni and iV; = 5 in the present 
work. We note that in SACRA and SACRA-ND, six grid 
points are prepared outside the refinement boundaries for 
buffer zone in each refinement level, following [4^] . The 
relation Li = 2L;+i is imposed for SACRA-ND. Irrespec- 
tive of we assign 2N + 1 vertex-centered grid points for 
[—Li : Li] where N was chosen to be 30, 40, and 50 for 
checking convergence: We monitored the violation of the 
Hamiltonian and momentum constraints, conservation of 
gravitational mass, and conservation of the area and spin 
of BHs (for stable model) , and checked that the conver- 
gence with improving the grid resolution is achieved. The 
grid spacing is Li/N and L^v, is chosen to be 1.2^^/^''"'^) 
irrespective of the number of d. For the best resolved 
run, the apparent horizon radius is covered by 30 - 40 
grid points for rapidly spinning BHs: We note that the 
initial coordinate radii of the BH horizon, r;i//i^/(''~^\ 
are in the range 0.584 - 0.630 (g = 1.143 - 0, minimum 
at q - 1) for d = 6, 0.750 - 0.707 (q ^ 1.013 - 0) for 
d = 7, and 0.812 - 0.758 {q = 0.832 - 0) for d = 8, re- 
spectively. Because of our choice of the puncture gauge, 
this coordinate radius increases by a factor of 1.3 - 1.5 
during evolution. 

A steep gradient always appears for geometric quanti- 
ties near the origin for rapidly spinning BHs in the punc- 
ture gauge, and this often causes a problem for stably 
evolving the BHs, although such region is not important 
for studying the stability of the BH at all. Thus, we em- 
ploy a very simple excision method for discarding this 
region. Specifically, for r < r^x, we set jij — Sij and 



A,j ^ K = 



Other quantities, a, /3 , x, and F*, are 
solved without any prescription. In this work, rex is cho- 
sen to be 0.3/x^/(''~^) or 0.4/i^/(''^''' which is typically 
^ 30 - 40% of the apparent horizon radius; thus, the 
excised region is well inside the BH horizon. Because no 
information escapes from the BH horizon and there are 
many grid points covering the inside of apparent horizon, 
the results are insensitive to the choice of Tex as far as 
it is sufficiently small. We note that this method is ac- 
ceptable for the case that an unstable BH deforms by a 
moderate degree (for rj < 0.5; see Sec. lIVBl for definition 
of rj). For a highly nonlinear deformation, however, we 
will have to develop other prescriptions for handling the 



7 



steep gradient near the central region, which is beyond 
the scope of this paper. 



B. Numerical results 

1. Critical spin for bar-mode instability 

During numerical simulation, we determine the appar- 
ent horizon and calculate its area and circumferential 
radii to investigate the properties of the unstable BHs. 
We also extract gravitational waves of quadrupole mode 
in a local wave zone. 

To determine the stability of a BH against bar-mode 
deformation, we monitor two quantities. One is a de- 
formation parameter of the BH horizon. To define this 
parameter, we first calculate circumferential radii of the 
apparent horizon along several meridians. Specifically, 
we measure the proper length of the meridians for (p = 
(and tt), 7r/4 (and 57r/4), tt/2 (and 3tt/2), and 37r/4 (and 
77r/4); the proper length of each meridian for a given 
value of ip{< tt) is defined by 



7r/2 



(^Vleei^p) + Vlseiv^ + ir^jdO. (41) 



Then, we define a bar-mode deformation parameter 



77 := 



Iq + Itt/2 



(42) 



The value of rj is zero for an axisymmetric apparent hori- 
zon, and increases as the deviation from the axial sym- 
metry is enhanced. Thus, rj is an indicator of nonax- 
isymmetric bar-mode deformation. A drawback of the 
diagnostic with this quantity is that the apparent hori- 
zon is a coordinate-dependent notion, and thus, rj may 
not exactly trace the deformation of the BH, although 
the apparent horizon agrees with the event horizon for 
a stationary spacetime, and hence, for 77 ^ 1, they are 
likely to agree at least approximately each other. 

The other more physical quantity for measuring the de- 
formation of a BH is gravitational waveform observed in a 
wave zone. This quantity should be coordinate-invariant, 
and thus, tells us whether the BH spacetime is stable or 
not with no ambiguity. In this work, we monitor a di- 
mensionless form of gravitational waves 



^yy 



-| {d~2)/2 



il/(rf-3) 
-| (d-2)/2 



l/(d-3) 



(43) 
(44) 



along the z-axis. Here, r is the coordinate distance from 
the center, and /i+ and /ix are regarded as the plus (+) 
and cross (x) modes of quadrupole gravitational waves 
in the wave zone. 



Figure [Tfa) plots the evolution of 77 for d = 6 and 
for the initial spin = ai/ y}/^ — 0.674 - 1.039 with 
A = 0.005. This shows that for i/^^^^ ^ 10, the value of 
77 increases or decreases in an exponential manner as 



tlT 

V = Voe' , 



(45) 



where r is the growth timescale and 7^0 is a constant. Fig- 
ure[IJb) plots the growth rate, , in units of /i^^/"^ as a 
function of g^. We determine r by the least-square fitting 
of the curve of In 77 with linear lines for 15 < tj [i^l'^ < 100 
or for the data of 15 < t/fj,^^^ and 77 < 0.1. These fig- 
ures show that for qi > 0.75, 77 grows exponentially with 
time, and otherwise, it damps exponentially. Therefore, 
the BHs with qi > 0.75 are dynamically unstable against 
bar- mode deformation [the critical value, gcrit , is found to 
be « 0.743 by interpolation; cf. Fig.[IJb)]. Because they 
are dynamically unstable, such BHs have a quasinormal 
mode with a negative value of the imaginary part of eigen 
angular frequency, w/ := Im(wQN), and the bar-mode de- 
formation grows in proportional to e'""'"' oc e~'^'* . On 
the other hand, the BH with qi < 0.74 is stable, and 
thus, uJi are positive for any mode (strictly speaking, 
there might exist a mode of w/ < with |a;/| less than 
^ 0.01/i~^/'^ because we performed the simulations only 
for a finite duration of order lOO^^/'^). 

Figure IHb) shows that the growth rate of the instabil- 
ity, r~^, increases monotonically with the spin g, and 
approximately linearly in q for 0.74 < q < 1. The 
growth rate (t^^ or ojj) is approximately described as 
Criq — qctit) where Cr is a constant « 0.51iJ,~^^^. 

The growth timescale of the bar-mode deformation 
for the BHs of spin slightly larger than gcrit is of order 
100/1^/"^ (the damping timescale for a stable BH with q < 
(7crit is also of this order). This is much longer than the 
spin period of the BHs defined by 27r/f7H = ©(lO/i^/^) 
where fl^L = a/(^+ + a^) is the angular velocity of the 
BH horizon. Thus, the instability grows slowly with 
a timescale much longer than the dynamical timescale 
^ r+/c for q ~ (?crit- For a much larger value of q, how- 
ever, the growth timescale is much shorter: We find that 
T is as short as 27r/riH for q ~ 0.95, and as Tr/fin for 
g 1.15. Simply extrapolating the approximately linear 
relation for r~^(g), r is likely to be much shorter than 
27r/riH for q ^ 1. Implications of these properties for 
the evolution of rapidly spinning BHs will be discussed 
in Sees. V B and VI. 

Figure [2l^ a) plots the longterm evolution of rj for d = 6 
and for = 0.801 with several magnitudes of initial per- 
turbation {A = 10-6, 0.001, 0.005, and 0.02) and for dif- 
ferent grid resolutions (A^ = 30, 40, and 50 for A = 0.005, 
and iV = 40 and 50 for A = 0.02). This shows that 
the growth rate is independent of the initial perturbation 
magnitude, implying that the instability sets in irrespec- 
tive of initial perturbation; even from an infinitesimally 
small perturbation, the instability grows spontaneously. 
The growth rate depends weakly on the grid resolution; 
for a poor grid resolution {N = 30), the growth rate 
is underestimated. However, for TV > 40, it appears to 



8 




(a) (b) 
FIG. 1: (a) Evolution of deformation parameter t] for d = 6 and for the initial spin qi = aj ~ 1.039, 0.986, 0.933, 0.878, 
0.821, 0.801, 0.781, 0.761, 0.750, 0.740, 0.718, and 0.674 (from the upper to lower curves) with A = 0.005. (b) The growth rate 
of r\, 1/t, in units of as a function of q (solid curve). The dashed curve denotes Q.^/2n. For qi > 0.75, the value of r] 

increases exponentially with time, and otherwise, an exponential damping is seen. For qi = 0.750 (thick solid curve in panel 
(a)) and 0.740 (below the curve of qi = 0.750), the growth and damping rates of r] are quite small, indicating that these BHs 
are close to the marginally stable state. 



approximately converge, and thus, the simulations with 
TV > 40 (grid spacing smaller than 0.03//^/"^) have an 
acceptable resolution. 

After the value of rj reaches the maximum, the growth 
of the instability terminates, and then, rj damps expo- 
nentially with time. The reason for this is that dur- 
ing the nonlinear growth of the bar-mode deformation, 
gravitational- wave emission is enhanced, and energy and 
angular momentum of the BHs are carried away. As a 
result, the BH spin parameter, decreases and at the 
saturation, it reaches a stable state with q < gcritj for 
which the sign of the imaginary part of eigen angular fre- 
quency for the corresponding quasinormal mode changes 
to be positive. The evolution of the BH spin will be ana- 
lyzed in detail in Sec. IIVB 21 Here, we show an indirect 
evidence that this interpretation is correct: The satu- 
ration value of 77 is larger for the BHs of larger initial 
spin, because the required spin-down fraction is qi — (/crit- 
Figure mb) plots the evolution of 77 for A — 0.02 and for 
relatively small initial spins, qi = 0.821, 0.801, and 0.781. 
Figure is the same as Fig. HJb) but for large initial 
spins qi = 0.878, 0.933, 0.986, and 1.039 with A = 0.005. 
These clearly illustrate that for the larger initial spin, 
a BH has to emit more gravitational waves for spinning 
down to reach a stable BH, and hence, the nonlinear 
growth of the bar-mode deformation has to continue un- 
til a high saturation value of 77 is reached. Figure [^^d) 
plots the maximum value of 77, 77niax, as a function of qi. 
This clearly shows that the value of 77max increases sys- 
tematically with qi'. The relation between 77„iax and qi is 
approximately written in the form (see Sec. VI B for 
an approximate derivation of this relation). This result 
suggests that for ^ 1, the BH may reach a highly de- 
formed state with 77 ^ 1 (although we were not able to 
evolve such ultra spinning BHs in the present work.) 

The numerical results obtained in this paper are ap- 



proximately derived using a semi-analytic calculation, as 
illustrated here. Based on the analytic calculations, we 
will show in Sec. |V] that the BH spin should be decreased 
by gravitational radiation reaction. 

It is worthy to note that for the larger initial spin qi ~ 
1, the damping timescale for the bar- mode deformation 
after the saturation is reached depends weakly on the 
initial spin: The timescale (the duration until 77 < 10"'^ 
is reached) is ^ IQQy}/^. This seems to be due to the 
fact that the damping timescale depends primarily on 
the BH state after the saturation is reached; i.e., the 
property of a BH of spin q < gcrit determines the damping 
timescale. The damping time duration is fairly long ~ 
100//^/"^, implying that the unstable BHs always have to 
spend a long time until they reach a stationary, stable 
state. 

Figure plots the evolution of rj for d ~ 7 and 
qi = = 0.719 - 0.960. This shows that for 

Qi ^ 0.735, rj grows exponentially with time, and oth- 
erwise, it damps exponentially. Figure [3Jb) displays the 
growth rate, t~^, in units of and shows that it 

monotonically increases, approximately linearly with q 
for 0.719 < 9 ^ 1 as in the case of d = 6. The critical 
spin for the onset of the bar-mode instability is deter- 
mined as (7crit ~ 0.730, and near q — q^it the growth 
rate (t~^ and uji) behaves as Cr{q — qmt) where Cr is a 
constant ~ 0.54/i~^/^. It is interesting to note that the 
coefficient Ct/x^^'''^""^-' for c? = 7 is close to that for d = 6. 

The critical spin, (/crit, for the onset of the bar- mode 
instability is also very close to that for d = 6. Simulations 
for d = 8 also clarified that the critical spin is Qcrit ~ 0.77, 
where q — aj for d = 8, again close to those for d — % 
and 7. This suggests that the critical value of q depends 
weakly on the value of d as long as d > 6 (note, however, 
that for d = 5, we found it a relatively large value, 0.87 
[H). The values of gcrit and the corresponding values of 
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FIG. 2: (a) Evolution of deformation parameter for d — 6 and for qi = 0.801 with initial perturbation amplitude A — 0.02, 
0.005, 0.001, and lO"'' (from the upper to lower curves) and with A'^ = 30 (dotted curve), 40 (dashed curves), and 50 (solid 
curves), respectively, (b) The same as (a) but for qi — 0.821, 0.801, and 0.780 (from left to right) with A = 0.02 and with 
A'^ = 40 (dashed curves) and 50 (solid curves), (c) The same as (b) but for qi — 1.039, 0.986, 0.933, and 0.878 (from the upper 
to lower curves) with A — 0.005. (d) The maximum value of 77 as a function of qi. The dashed line and solid curve denote 
??max ~ ^{qi — Qcrit) aud rclatiou (|68p . respectively (see Sec. VI B for an approximate derivation of these relations). 



TABLE I: The values of Qcrit and the corresponding values of 
a, and Cp/Ce for d = 5 - 8. Cp/Ce is defined in Eq. (|48t . 



d 


5 


6 


7 


8 


9crit 


0.87 


0.74 


0.73 


0.77 


a* 


1.76 


0.91 


0.83 


0.86 




0.38 


0.65 


0.68 


0.67 



a* = are summarized in Table HI 

As in the case d = 6, the growth rate, t^^, increases 
monotonically with the value of qi for d = 7. For qi > 
0.92, it is larger than the BH spin frequency H.u/'^tt. This 
indicates the growth timescale is much shorter than the 
BH spin for ^ 1 irrespective of the dimensionality. 

Figure ^c) plots the longterm evolution of 77 for d = 7 
with higher spins, qi = 0.813 - 1.017. For these BHs, 
the growth of the bar-mode deformation saturates at a 
large value of r] as r?max = 0.1-0.5 due to gravitational 
radiation reaction, and then, rj decreases exponentially 
as in the case of d = 6. At the saturation, the BH spin 



parameter, q, is likely to be slightly smaller than gcrit- As 
in the case of d = 6, the damping rate of the bar-mode 
deformation after the saturation is reached depends only 
weakly on the initial value of qi. The reason for this is 
that it depends primarily on the spin achieved just after 
the saturation is reached {q < qait), as mentioned before. 
One point to be noted is that the damping timescale is 
much longer than the growth timescale for qi ^ 1 (in the 
6D case, two timescales are not very different and smaller 
than 50 fj}^^). The damping timescale is > lOOn^^'^ in the 
7D case. This suggests that it always takes a long time 
for the unstable BHs to reach a stable state for d> 7. 

To confirm that the bar-mode deformation of the BHs 
is indeed physical (the growth of rj is not due to a coor- 
dinate choice for finding apparent horizons), we extract 
gravitational waves in the wave zone. Figure|3]plots grav- 
itational waveforms, and h^, as a function of a re- 
tarded time t — r ioi d ~ Q and for qi « 0.801 and 0.986. 
We note that gravitational waveforms for different values 
of qi are qualitatively similar. Here, r denotes the coordi- 
nate distance of extracting gravitational waves. We also 
plot together evolution of 77/2 as a function of t. As in 
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FIG. 3: (a) The same as Fig.[T]but for d = 7 and for qi = a/^i^/* = 0.960, 0.903, 0.844, 0.813, 0.783, 0.767, 0.751, 0.735, and 
0.719 (from the upper to lower cmves) with A'^ = 50. For qi > 0.73, the value of i] increases exponentially with time, otherwise, 
an exponential damping is seen. For q = 0.735, the growth (or damping) rate of i] is close to zero, implying that this BH is 
close to the marginally stable one. (b) The growth rate of r], 1/r, in units of /i-i/^ as a function of q (solid curve). The dashed 
curve denotes JIh/Stt. (c) The same as panel (a) but for long runs with qi — 1.017, 0.960, 0.903, 0.873, 0.844, and 0.813. 



the behavior of 77, the amphtude of gravitational waves 
increases exponentially with time in the early phase, and 
after the saturation point is reached, it starts exponential 
damping. For the larger initial spin, the growth timescale 
is shorter; for qi = 0.986, the grow timescale is compa- 
rable to the oscillation period. h+ and behave essen- 
tially in the same manner except for a phase difference 
of 7r/2. The amplitude of gravitational waves is approxi- 
mately equal to 77/2 for the small amplitude and slightly 
smaller than 77/2 for the amplitude {h^ + hl,y/^ > 0.1. 
This slight disagreement in the amplitude for the high- 
amplitude case is likely due to the fact that the nonlinear 
deformation of the BH enhances the amplitude of gravi- 
tational waves of modes other than the quadrupole ones, 
suppressing the quadrupole modes. It is worthy to point 
out that the oscillation frequency remains approximately 
constant (besides a small secular shift associated with the 
evolution of the BHs). This indicates that one fundamen- 
tal quasinormal mode contributes to the instability. 



2. Evolution of spin and area 

Now, we illustrate how the BH spin evolves as a result 
of gravitational radiation reaction. To approximately de- 
termine the value of the BH spin from the shape of appar- 
ent horizon, we calculate the ratio of a polar circumferen- 
tial length, Cp, to the equatorial circumferential length, 
Ce, for which we define 

- (46) 

Ce / /r^rfy', (47) 
Jo 

where the integral for Ce is performed along the surface 
of the apparent horizon at = tt/2. For the Myers- Perry 
BH of single spin, the ratio Cp/Ce is given by 

which is a monotonically decreasing function of q = 
u/^i/('i-3) fQj. g^jjy given number of d (see Fig.[S]). Thus, 
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FIG. 4: (a) and (b): + and x modes of gravitational waveform (solid curve) from an unstable BH for d — 6 and for qi = 0.801 
as a function of a retarded time defined hy t — r where r is the coordinate distance from the center. We also plot rj/2 as a 
function of t (dashed curve). The initial condition is A = 0.005 and the result for the grid resolution of A/" = 50 is plotted, (c) 
and (d): The same as (a) and (b) but for d = 6 and for qi = 0.986 with A = 0.005. 
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FIG. 5: Cp/Ce as a function of g = a/ fi^^^'''^'' for d = 5 - 8. 



as far as the state of a BH is close to an axisymmetric 
stationary BH, we may use it for measuring the spin, q. 
We note that 5D BHs become dynamically unstable only 
for a small value of Cp/Ce < 0.38 d, whereas for d > 6, 
the critical values of Cp/Ce are rather large universally; 
Cp/Ce = 0.65 for d = 6, 0.68 for d=7, and 0.67 for d = 8 



(listed in Table I). This illustrates that the 5D BHs are 
qualitatively different from the higher-dimensional BHs 
with d > 6. 

Figure [SJa) plots the evolution of Cp/Ce for d = 6 
and for relatively small initial spins qi ~ 0.821, 0.801, 
and 0.781. The initial value of this ratio is w 0.587 for 
q, = 0.821, sa 0.602 for q, = 0.801, and w 0.618 for 
qi — 0.781, respectively. When rj is much smaller than 
0.1, the value of Cp/Ce remains approximately constant. 
With the nonlinear growth of the bar-mode deformation, 
it starts increasing due to the decrease of the spin by 
gravitational radiation reaction, and it eventually set- 
tles to a constant « 0.700 for qi = 0.821, « 0.690 for 
qi — 0.801, and w 0.677 for qi = 0.781, respectively, 
when 77 <C 1 being achieved (we adopt the values of 
Cp/Ce when we stopped the simulations). We note that 
with a poor grid resolution, the value of Cp/Ce spuri- 
ously increases even for 77 <g; 1 because numerical dis- 
sipation tends to decrease the BH spin. This spurious 
variation is suppressed for the better grid resolutions (see 
also Fig. [7]). The final values of Cp/Ce indicate the final 
spin qf sa 0.675, 0.688, and 0.705 for q, = 0.821, 0.801, 
and 0.781, respectively, for which the BHs appear to be 
stable against any perturbation. It is worthy to note that 
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FIG. 6: (a) Evolution of Cp/Ce for d = 6 and for relatively small initial spins qt = a/fi^'^ = 0.821 (Cp/Ce ~ 0.587 at t = 0), 
0.801 {Cp/Ce « 0.602 at t = 0), and 0.781 {Cp/Ce « 0.618 at t = 0) as a function of time. For qt = 0.801, the results with 
A = 0.02 and 0.005 are plotted. For qi — 0.821 and 0.781, the results with A = 0.02 are plotted. The solid and dashed curves 
denote the results for A'^ = 50 and 40, respectively. The thin dotted line denote Cp/Ce = 0.647 which is the value of Cp/Ce for 
q = gcrit. For qi = 0.821, we stopped the simulation at t/fi^^^ ~ 370 because the BH reaches approximately stationary state, 
(b) The same as (a) but for the large initial spins g,; = 0.878, 0.933, 0.986, and 1.039 with A = 0.005. Cp/Ce « 0.542, 0.499, 
0.460, and 0.422 at t = 0, respectively. 
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FIG. 7: The same as Fig. [6] but for the area of apparent 
horizon for qi — 0.801 with A'^ = 30 (dotted curve), 40 (dashed 
curves), and 50 (solid curves) and for qi — 0.781 and 0.821 
with A'' = 50 (lower and upper dot-dashed curves). Ao is the 
initial value of the area, Qd-2fJ'r+. 



the larger initial spin leads to the smaller final spin be- 
cause a larger amount of gravitational-wave emission is 
enhanced during the evolution (see also Sec. |V] for the 
reason to this); the value of qf is approximately written 
as qj = 2gcrit — 9i for qt = 0.781 - 0.821. 

Figure [6ljb) plots the evolution of Cp/Ce for d = 6 
and for larger initial spins qt = 0.878, 0.933, 0.986, and 
1.039. As in the case of relatively small initial spins, the 
value of Cp/Ce increases with the evolution, and even- 
tually settles to a constant after the BH is stabilized 
by gravitational radiation reaction. However, the final 
value of Cp/Ce depends only weakly on the initial spin, 



and hence, the final spin estimated from Cp/Ce is much 
larger than the value of qj = 2qcrit — qi. Figure[Sl[b) shows 
that the fractional change of Cp/Ce after the BH is sta- 
bilized is smaller than that before the stabilization. This 
indicates that the fraction of spin-down is suppressed af- 
ter the evolution of the nonlinear bar-mode deformation 
is saturated; again, this seems to be due to the fact that 
the evolution of the BHs after the saturation is achieved 
is determined by the property of the BH of spin slightly 
smaller than gciit- For qi = 0.878 - 1.039, the final spin 
is relatively in a narrow range as g/ ~ 0.65 - 0.61. The 
results we found suggest that even for qi ^ 1, the final 
spin qf will not be close to zero, but moderately large as 
~ 0.6. 

Qualitatively the same results for the final spin are 
obtained for d ~ 7. In this case, Cp/Ce relaxes to 0.72 ± 
0.01 for 0.85 < qi < 1, implying that the final spin is in a 
narrow range qf = 0.66 - 0.69; even for qi S> 1, the final 
spin qf will not be close to zero, but moderately large as 
~ 0.65. The results of qf for d — 6 and 7 together with 
the weak dependence of gcrit on d suggest that also for 
d > 8, the final spin will be moderately large > 0.6 - 0.7. 

The area of the BH has to increase as a result of spin- 
down [l^l • Figure [7] plots the evolution of the horizon 
area for d = 6 and q^ = 0.781, 0.801, and 0.821 with 
N = 50. For qi — 0.801, we plot the results for differ- 
ent grid resolutions. Figure [7] shows that the area indeed 
increases. As noted above, the area gradually increases 
with time even in the state of 77 1. This is partly due 
to numerical dissipation. Indeed, the error decreases sig- 
nificantly with improving the grid resolution; for A'^ = 50, 
the spurious increase of the area seems to be negligible, 
and hence, we can determine the magnitude of the physi- 
cal increase due to gravitational-wave emission. The area 
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increases by « 1.0, 1.4, and 1.8% for qi = 0.781, 0.801, 
and 0.821, respectively. For the larger initial spins, the 
area increases more due to gravitational radiation reac- 
tion. A noteworthy fact is that the fractional increase 
of the area is much small than those for spin and mass 
energy of the BH (cf. Fig. H). 



V. PROPERTIES OF UNSTABLE BLACK 
HOLES 

A. Conditions for spontaneous gravitational wave 
emission 

The first law of BH thermodynamics allows us to de- 
termine the variation in the BH area SA as 13311 



87rGd 



(49) 



where k is the surface gravity of the BH horizon, k = 
[2r'l^^ + (d- 5)fi\/{2fir+), which is positive for the BHs 
considered in this paper. SE and SJ are the variations 
of energy and angular momentum of the BH, and if the 
energy and angular momentum are carried away by grav- 
itational waves, they should be negative. If we assume 
that gravitational waves of monochromatic wavelength 
are emitted, the following relation holds; 



SE = —SJ < 0. 
m 



(50) 



Here, ui is the (real) angular frequency of the unstable 
mode [we regard it as w = Re(a;QN) in the following], 
and m is the azimuthal quantum number for which we 
set m = 2 because we focus on the bar-mode instabil- 
ity. The assumption of monochromatic wave emission is 
approximately correct because gravitational waves asso- 
ciated with the fundamental quasinormal mode are likely 
to be most strongly emitted. Indeed, Fig.|4]indicates that 
this is the case. 

Substituting the relation ([SU]) into Eq. P^t . we obtain 



needs an artificial wave injection which satisfies the con- 
dition A: < 0. We note that for the spontaneous emission 
of gravitational waves, one further needs to require the 
presence of a quasinormal mode with negative imaginary 
part for its angular frequency, w/ < 0. 

To confirm that the superradiance condition is satisfied 
for the oscillation mode of unstable BHs, we first deter- 
mine the real value of cj by performing the Fourier trans- 
formation of /i+ for the BHs. We identify the angular 
frequency at the spectrum peak as cj. Figure [5] plots the 
resulting values for uj/2 (points) as well as i^H (curves) 
as functions of q. We also plot the results for 5D BHs ob- 
tained in the previous paper [l|. This figure shows that 
the unstable modes with w/ < always satisfy the su- 
perradiance condition flu > uj/2 and gravitational waves 
can be spontaneously emitted from the unstable BHs. A 
noteworthy point is that W(U^/('^^'^''(« 1.0-1.1) depends 
weakly on the spin g; for d = 6, uj decreases slowly with 
q for 0.7 < q ^ 1-1 and for d — 7 and 8, u behaves in a 
manner similar to that of d = 6 for 0.7 < g < 1. 

In the following, we refer to the value of q, at which 
the real part of the quasinormal frequency becomes equal 
to the threshold frequency for the superradiance, as qsR- 
The superradiance condition is satisfied for the quasi- 
normal modes when q > qsR, where ggR ~ 0-6 for d = 5, 
gsR ~ 0.65 for d = 6, qsR ~ 0.7 for d = 7, and gsR ~ 0.75 
for d = 8. Because the maximum value of V.}ifj,^/^'^~^^ de- 
creases as the number of d is increased and w/z^^'-''^"^-' for 
a given value of q depends weakly on d, the value of ggR 
will increase as d is increased. This indicates that the 
critical value of spin, ^crit , for the onset of the bar- mode 
instability (i.e., for the spontaneous gravitational- wave 
emission) will also increase with d for d > 8. 

We note that Fig. [8] shows fin ~ Lu/m <^ iln for 6 < 
d < 8, and thus, the area increases only slightly even if 
a large amount of angular momentum is dissipated by 
gravitational- wave emission [see Eq. ([ST]) ]. This explains 
why the area increases only slightly (see Fig. [7]) even when 
q changes by a large factor ^ 0.1 - 0.2. 



^SA=(nu--)\sj\. 



(51) 



This shows that only for Oh > uj/m, 5 A becomes positive 
and the evolution by emission of gravitational waves is 
allowed without violating the area theorem [20|. This is 
the so-called superradiance condition [i^ : If this condi- 
tion holds, the energy fiux of ingoing waves of an angular 
frequency k = lo — mfin a-t the BH horizon becomes neg- 
ative because A: < 0. 

This is a necessary condition that the unstable quasi- 
normal mode should satisfy for the spontaneous emission 
of gravitational waves. However, this is not the sufficient 
condition. The superradiance is similar to secular in- 
stability in the terminology of an instability for rotating 
stars [U [13] , because for the superradiance emission of 
gravitational waves from a stable BH with luj > 0, one 



B. Mechanism of spin-down 

Next, we consider the evolution of nondimensional spin 
parameter, q, as a result of gravitational-wave emission. 
The variation for this quantity, in the assumption of 
monochromatic gravitational- wave emission, is written as 

Sq ^ 8nG,f,-^''-^y(''-'^n,l,(^SJ - ^SE^ 

= ^8nG,f,-('-'y('-^^n,l,\6J\ (l - (52) 



^ \SJ\ ( 2 aLo 
5q = -j-q 1 - — r — 

J \ d — 6 m 



(53) 
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FIG. 8: as a function oi q = a/^^^^'''~^^ for d = 5 - 8 (curves), and uj/m = iu/2 for selected values of the spin parameter 
for d = 5 - 7 (points). Tiie units of the vertical axis are /i"^''''*"^' . We omit the points for d = 8 because they are located in 
the vicinity of the points for d = 6 and 7: The values of ujfi^^^''~^^ /2 are in a narrow range 0.49 - 0.55 for the chosen values 
of q and for d = 6 - 8. The squares on the curve of Qh approximately denote the critical points for the onset of bar-mode 
instability (i.e., for spontaneous gravitational-wave emission). 
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FIG. 9: \Sq\ (left) and 5E/M (right) as functions of t for d = 6. 
qi = 0.781 and 0.821 with A = 0.02 are displayed for iV = 50. 
because the BH reaches approximately stationary state. 
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The results for qi = 0.801 with A = 0.005 and 0.02, and for 
For qi = 0.821, we stopped the simulation at t/fj}^^ « 370 



As found from Eq. (|51l) , the gravitational- wave emission 
is possible only for flu > uj/m. Thus, 



2 auj 
1 - - — - — > 1 



d — 3 rn 



d-3 



1 



rf- 3: 



(54) 



Because 2/(d - 3) < 1 and aV(?-+ + a^) < 1 for any BH 
with d > 5, we find that i5g is always negative for the su- 
perradiance mode. Thus, the unstable BHs have to spin 
down by gravitational radiation reaction and evolve to- 



ward a stable state, as we found in numerical simulation 
(cf. Fig.©. 

We quantitatively confirm from Eq. ([55)) that the spin- 
down of the unstable BHs is indeed due to gravitational- 
wave emission as shown in the following. Using a for- 
mula of gravitational- wave luminosity [51] (the Landau- 
Lifshitz pseudo tensor [i^l also gives the same formula 
after partial integration), the luminosity of quadrupole 
gravitational waves {I = \m\ = 2 mode) is written as 



dE 



WnGd 



(d-2)(d+l) 



(d-2)/(d-3) 



.(55) 
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For monochromatic gravitational waves with angular fre- 
quency uj, the luminosity may be rewritten as 



dE_ _ iu^{hl + hl)nd-2 



{d - 3)rf 



id-2){d+l) 



(d-2)/(d-3) 



(56) 



In the following, we assume /i+ [uj{t — r)] = hx [(^{t — r) + 
tt/2] (and set h = ^Jh\ + h\), and w =const for simplic- 
ity, which approximately hold as mentioned before. 
Rewriting Eq. (|52p in the form of spin evolution as 



dq 
'dt 



B7rmGdAi-('^-2)/(d-3) 



d-2 



dE 



dt 



d ~ 3 m 



,(57) 



and integrating this equation, we may infer the spin-down 
history of a BH for the whole evolution process by 



1 c- / M , . mujh? d(d — 3) 



1 - 



d — 3 m 



(58) 



We also calculate dissipation history of the mass energy 
by 



SEit)^ldt'§, 



(59) 



where we use Eq. ((56| for dE/dt. We note that Eq. ([52)) 
is valid only for the case that deformation from a station- 
ary BH solution is small. Thus, we here apply Eq. (|58|) 
only for the small values of qi — Qcrit < 0.1 for which 
the maximum value of is < 0.15 and the deformation 
remains in a weakly nonlinear level. 

Figure [9] plots \5q\ and 5E/M as functions of t for 
d = 6 and for q^ = 0.821, 0.801, and 0.781. Here, we 
assume that auj = 1.05gi {ujjj}/'^ = 1.05) taking into ac- 
count the results shown in Fig. HI and also neglect the 
dissipation of mass and spin of the BHs in computing 
the integrals of Eqs. (|58l) and (|59l) . because the variation 
of these quantities gives a minor effect on 6q and SE for 
Qi — QcTit < 0.1; the error is likely to be less than 10%. 
As shown in Fig. [HI we obtain (Sgtot ~ 0.14, 0.11, and 
0.07 for q, w 0.821, 0.801, and 0.780, respectively, which 
agree with the results determined from the final value of 
Cp/Ce within the error A{5q) « 0.006. Based on this 
good agreement, we conclude that the gravitational radi- 
ation reaction determines the spin-down of the unstable 
BHs and that in our simulations, the spin-down process 
is computed accurately. We also note that the numerical 
results for the fractional increase of the area (cf. Fig. [7]) 
agree well with that calculated from 5q and 5E obtained 
here. 

As we have described above, the evolution process of an 
unstable BH with g < 1 is determined by gravitational ra- 
diation reaction. As far as the growth rate of the unstable 
mode is comparable with or smaller than the spin-down 
rate due to gravitational radiation reaction, the scenario 
presented here will be always correct. Figures [IJb) and 
[3]Jb) show that for the unstable BHs with qi < 1.1, which 



have the growth timescale (r) longer than Tr/fin, this sce- 
nario holds. We, however, note that this scenario might 
not be valid ior q ^ 1: If the growth timescale is too 
short (e.g., toj < 1), the BH could not emit a significant 
amount of gravitational waves in t for spinning down, 
and as a result, it can achieve a state for which the de- 
formation parameter is of order unity. Extrapolation of 
the results shown in Fig. [UJd) indeed suggests that for 
Qi }^ 1-6, 1] may reach > 1 for d = 6. For such case, grav- 
itational radiation reaction may not prevent the growth 
of the bar-mode deformation, and then, highly nonlinear 
deformation may lead to fragmentation of the BH (in the 
classical argument, formation of a naked singularity may 
be the result), as discussed in [s^l- In the present work, 
we have not pursued the possibility that rj reaches ~ 1, 
because a longterm simulation for such an ultra spinning 
BH is not an easy task technically. We leave such a work 
for the future. In Sec. IVI B| we speculate a condition for 
an ultra spinning BH to cause fragmentation. 

Before closing this section, we note that the fraction 
of total radiated energy of gravitational waves is much 
smaller than that of angular momentum. The reason for 
this is that the following relation holds: 



5J_d-26E 
Jq au! Mq 



(60) 



Here, Mq and Jq are initial mass and angular momentum 
of a BH. For d = 6 and qi = 0.801, for example, SE/Mq 
is only w 4%, whereas SJ/Jq is ~ 20%. This significant 
angular-momentum dissipation is essential for the spin- 
down of the unstable BHs. It is worthy to note that 
the angular momentum dissipation rate is even larger for 
the larger number of d for a given value of q because 



LUfJ, 



l/(d-3) 



depends weakly on d for the unstable BHs. 



VI. SUMMARY AND DISCUSSIONS 
A. Summary 

We show by numerical-relativity simulation that BHs 
spinning sufficiently rapidly are unstable against nonax- 
isymmetric bar-mode deformation for d = 6, 7 and 8, as 
in the case d = 5 [l[. In this instability, gravitational 
waves are emitted spontaneously. The critical BH spin 
for the onset of the bar-mode instability is ^crit ~ 0.74 
for d = 6, « 0.73 for d = 7, and « 0.77 for d = 8, re- 
spectively. Thus, the critical value is smaller than unity, 
and depends only weakly on the dimensionality for d > 6. 
After the instability sets in, the degree of the bar-mode 
deformation increases exponentially with time, and even- 
tually, it saturates. The saturation is caused by gravita- 
tional radiation reaction. For the larger initial spin, the 
growth timescale is shorter. The degree of the maximum 
deformation is larger for the larger initial value, approx- 
imately in proportional to qi — Qcrit- 

The unstable BHs emit gravitational waves signifi- 
cantly, and then, spin down, settling to a stable BH with 
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q < 9crit- The final value of the BH spin, g/(< qcHt), 
is smaller for the larger initial spin: For d — 6 with 
qi ~ 9crit < 0.1, the final spin is approximately written as 
qf = 2gcrit - qi, and for 0.1 < qi - ^crit < 0.3, the final 
spin is in a narrow range qf ~ 0.6 - 0.65. For d = 7, 
qf « 0.66 -0.69 for 0.85 < gi < 1. The smaller value of 
qf for the larger initial spin qi is due to a larger amount 
of gravitational- wave emission induced by the higher de- 
gree of the bar-mode deformation. The results for q.^ ^ 1 
indicate that for ^ 1, the final spin is likely to be 
^ 0.6 for d = 6 and ^ 0.65 for d = 7, if the BH settles to 
a stable state. 

The dynamically unstable BHs satisfy two conditions. 
One is that a quasinormal mode satisfies a superradiance 
condition; the real part of the eigen angular frequency 
for the quasinormal mode, uj = Re(a;QN), has to satisfy 
the condition uj < jnQn . The other is that the imaginary 
part of the eigen angular frequency, uji = Im(wQN), is 
negative. The superradiance condition referred to here is 
satisfied not only for the dynamically unstable BHs but 
also for the stable BHs with the spin smaller than the crit- 
ical value (gcrit)- The minimum spin for which the real 
quasinormal frequency satisfies the superradiance condi- 
tion, qsR, is ~ 0.6 for d = 5, ~ 0.65 for d = 6, ~' 0.7 
for d = 7, and ^ 0.75 for d = 8, respectively. Thus, gsR 
increases with the number of d. This suggests that qcnt 
also increases with d for d > 8, because gcrit has to be 
larger than qgji for the spontaneous gravitational-wave 
emission. 

It is worthy to note that the dynamically unstable BHs 
we found always satisfy the above two conditions. If the 
superradiance condition were not satisfied, a BH could 
not emit gravitational waves spontaneously, even if it is 
dynamically unstable (w/ < 0). Probably, there would 
exist a mathematical proof for the relation qcnt ^ 'ZsR - 

The value of gcrit depends weakly on the dimensionality 
for d > 6 as mentioned above. By contrast, this value is 
slightly larger for d = 5 (gcrit ~ 0.87) [l| and absent 
for d = 4. This shows that properties of the spinning 
BHs are qualitatively similar for d > 6, whereas those for 
d = 4 and 5 have their intrinsic properties. 



B. Axisymmetric vs nonaxisymmetric vs 
fragmentation instabilities 

The critical spin for the onset of the bar-mode insta- 
bility found in this paper is much smaller than that for 
the onset of axisymmetric instabilities [l^ll^, q > 1.56 
- 1.80 for d = 6 - 9, respectively (the minimum spins 
depend weakly on the dimensionality). Considering the 
analogy with the instabilities on rotating stars, wc find 
that this is a quite reasonable consequence [47, 52-54]: 
Rapidly and rigidly rotating stars can be unstable against 
axisymmetric instabilities (e.g., against toroid or ring for- 
mation) only for an extreme case in which the axial ratio 
of the polar axial length to the equatorial one is quite 
small (e.g., smaller than 0.171 for the incompressible fluid 



[52|). By contrast, a variety of nonaxisymmetric insta- 
bilities can set in even for the case that such ratio is not 
very small (for < 0.58, secular instabilities set in, and for 
< 0.31, dynamical instabilities set in for the incompress- 
ible fluid [13, [H]). Also, the bar- mode instability occurs 
at the lowest critical spin among many other instabilities 
for most rotating stars; the bar-deformation is the most 
efficient way for decreasing total energy of the system. 
The bar-mode instability may be most relevant for any 
self-gravitating spheroidal objects spinning rapidly. 

According to an estimate based on the BH thermo- 
dynamics by Emparan and Myers [ssj . the critical spin 
for the onset of a nonaxisymmetric instability is approx- 
imately unity, gfrag ~ 1, irrespective of the number of d 
for d > 6. The assumption in their argument is that a 
rapidly spinning BH of area Ai will fragment into two 
non-spinning boosted BHs of the total area A f which is 
larger than Ai. For the instability found in this paper 
for 9i ^ Ij the growth of the perturbation saturates at a 
weakly nonlinear level at 77 < 0.5, because gravitational- 
wave emission suppresses the further nonlinear growth. 
Thus, the fragmentation docs not occur for ^ 1. How- 
ever, this does not imply that the fragmentation instabil- 
ity is not relevant for any spin. As we showed in Figs. [T] 
and[31 the growth timescale of the unstable mode depends 
strongly on the magnitude of the spin, and for g > 1 , it 
is shorter than a spin period of the BH, r < 27r/r2H- For 
g ^ 1, the growth timescale is likely to be even shorter. 
In such case, gravitational-wave emission will play a mi- 
nor role and nonlinear growth will continue until a state 
of 77 ^ 1 is achieved, as discussed in Sec. V. Then, highly 
nonlinear bar-mode deformation may lead to a fragmen- 
tation of the BH due to the Gregory-Laflamme instability 
[3^ . In the following, we infer how large initial spin is 
required for achieving the fragmentation focusing on the 
case d = 6. 

As we showed in Figs. [1] and O the degree of bar- mode 
deformation for the unstable BHs increases exponentially 
with time as far as 77 < 0.1 as 

V = -, (61) 
r 

where r is a function of g, and for d = 6, sa 
Ct('? — <Zcrit) with Cr ~ 0.51^"^/'^. Gravitational ra- 
diation reaction plays a crucial role for suppressing the 
growth of the deformation. A nonlinearity associated 
with this effect plays an important role for larger values 
of 77. We take into account this nonlinear effect by phe- 
nomcnologically replacing t to r(l -|- Csri/2) where Cg 
is a constant of order unity. Because the amplitude of 
gravitational waves, h, is approximately written as 77/2, 
we may assume the following approximate equation for 
the growth of gravitational-wave amplitude: 



17 



Combining this relation with Eq. (j57|) . we obtain 



{1 + Csh) — 
aq 



4 (d-2)(d + l) 
Tmuj (d — 3)d 



1 



auj 



d — 3 m 



This equation is valid only for h <^ 1, but in the following, 
we use it even for h = 0(1) to infer the maximum value 
of ft, as a function of qi. For d = 6 with m = 2 and 
9crit < g < 1.1, this relation is approximately written as 



Because the factor, 1 — auj/3, is a slowly varying function 
of q in the range 0.76 - 0.64 for q = 0.7 - 1.1, we assume 
that it is a constant and denote it by F~^. In addition, we 
assume uj^^^^ w 1 for simplicity. Using the approximate 
relation for Cr, we finally obtain 

dh'^ 

(l+C,/l)— «-1.6F((J-(Jerit). 

dq 



(65) 



Here, we assume that the maximum value of h (/ipcak) 
is achieved when the spin q reaches gcrit as a result of 
spin-down. Then, integrating Eq. (j65p from q = q.^ to 
q = gcrit gives 



^pcak + -^^poak ~ 0.8F((7i - 9crit)^, 



or, using /ipoak ~ '7max/2, 



3.2F{q, - (Zcrit)' 



(66) 



(67) 



Here, 3.2F is « 4.2 - 5.0 for g « 0.7 - 1.1. 

Equation (|67p is indeed a good approximate relation 
for an appropriate choice of Cg . The solid curve of Fig. [5] 
(d) is 



7?2 
/max 



2??,nax = 4(9, - gf„it)' 



(68) 



and fits the numerical results well. It is found that for 
a high spin qi > 0.9, the nonlinearity associated with 
gravitational radiation reaction [the second term in the 
left-hand side of Eq. plays an important role in de- 
termining the value of jymax. 

Assuming that Eq. (|68| is valid even for > 1.1, we 
expect the maximum value of rj and find that r^max 1 
for qi > 1.6 (note that the possible maximum value of i] 
is 2). This suggests that the fragmentation may occur for 
an ultra spinning BH and the critical spin for the frag- 
mentation, (7frag, is ^ 1.6; this predicted critical value is 
much larger than the value predicted in [33.] . The critical 
spin for the onset of axisymmetric instabilities is fairly 
close to gfrag as gaxis ~ 1.56 for d = 6 [1^. Therefore, 
for the ultra spinning case, the BH may fragment in a 



complicated manner. We note that for r]^ 



> 1, 



TABLE II: The values of maximal impact parameter 
&max//i^^'''~'^' for the apparent horizon formation obtained in 
[5^ and the corresponding nondimensional spin of the system. 



d 


5 


6 


7 


8 


t-'max / 


1.24 


1.47 


1.59 


1.66 


^max 


0.93 


1.47 


1.98 


2.50 



role for the evolution of the deformed BH. For exam- 
ple, emission of gravitational waves with modes other 
than the lowest-order quadrupole mode could contribute 
to the spin-down and to suppressing the growth of the 
bar-mode instability. Thus, the critical spin for the frag- 
mentation may be even larger than 1.6. Finally, we note 
that essentially the same argument holds also for d> 7. 



C. Implications for mini black hole evolution 

The bar-mode instability found in this paper changes 
the hypothetical picture for the evolution of a mini BH 
which may be formed in particle accelerators. It is nat- 
ural to expect that most of the mini BH is formed for 
a large impact parameter, which is close to a critical 
value, &max, in two-particle collision, and that the BH 
should be rapidly spinning at its formation. In fact, 
analyses [ssl |56j for high- velocity two-particle collisions 
indicate that a BH can be formed for impact parame- 
ters such that the resulting BH is rapidly spinning with 
q > I for d > 6. To clarify this point, we first review the 
results of an analysis for high-velocity particle collision 
with the impact parameter b and with the energy of each 
incoming particle p [ssl l56j . For this phenomenon, the 
total mass energy and angular momentum of the system 
is M = 2p and J = bp, respectively, assuming that the 
particles move with the speed of light. As in the case 
of the Myers-Perry BH of single spin, we define spin and 
mass parameters, a and fi, for this system using the same 
formulas as Eqs. (|4]) and ([5]), and then calculate nondi- 
mensional spin q = a//i^/(''^'^). Table IhI shows the max- 
imal impact parameters ^max/M"'"/*-'^"'^-* for the apparent 
horizon formation obtained in [56j and the corresponding 
value of the spin q ~ (/max for rf = 5 - 8. The value of 
9max is larger than unity for d = 6 - 8, and thus, the 
resulting BH is likely to be rapidly spinning so that it 
can be subject to the dynamical bar- mode instability ^. 



eral additional nonlinear effects will play an important 



^ In this analysis, we assume that the mass and angular momentum 
for the formed BH are equal to those of the system for simplicity. 
However, this is not likely to be the case because gravitational 
waves should be significantly emitted during the collision. The 
value of Qmax is nothing but a roughly approximated value. As we 
showed in this paper, angular momentum will be more efficiently 
emitted by gravitational waves than energy. Thus, the resulting 
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The phenomenology of a mini BH formed in particle 
accelerators is determined by two timescales: One is the 
timescale of gravitational- wave emission tgw l/l'^/l- 
The other is the timescale of Hawking radiation th- For 
convenience, we denote these timescales in terms of the 
Planck mass mp . Among several manners of defining the 
Planck mass (summarized in [5|), we adopt the definition 
as 



mp 



(69) 



Then, the Planck time is defined by rp := l/mp in the 
natural units. 

The timescale of gravitational- wave emission is written 

as 



Tgw — (^Gw^p — 

TOP 



l/(d-3) 



(70) 



where Mq is the initial BH mass. If the formed BH is sta- 
ble, Tgw should be defined by l/|a;/|, and the coefficient 



<-^nw — — 



TT{d - 2)f7rf_2M 



l/(d-3) 



(71) 



For the 6D head-on collision {q = 0), for example, 
Cgw ~ 3.8 [i^], and this factor depends weakly on the 
dimensionality as long as g = 0. For spinning BHs, it 
becomes larger, and Cgw ^ 1 for q ~ gcrit as found in 
the present paper; practically Cqw = O(IOO), if gravita- 
tional radiation reaction is taken into account. 

If the value of q is further increased to be 9 > gcrit, the 
value of Cgw should be approximately written as 



Cgw = 



2n 



-I l/(d-3) 



7r(d - 2)fld-2lJ' 



C, 



GWO, 



(72) 



where the first term in the right-hand side is associated 
with the growth timescale for the bar-mode deformation 
and the second one is the timescale of damping after the 
saturation is achieved, which is approximately equal to 
Cgw in Eq. ([7T|) for q < gcrit and thus of order 10 - 10^ 
[cf. Fig. mb) and (c)]. The first term is comparable to 
the second term for gcrit < 9 ^ 9crit -1-0.1 but the second 
term dominates for the larger value of q. 
The timescale of Hawking radiation is 



m = Chtp 



Mo 
mp 



(d-l)/(rf-3) 



(73) 



where Ch is a ^-dependent constant and determined tak- 
ing into account all emission processes (i.e., emission of 
scalar, spinor, vector particles and gravitons). For the 



spin could be smaller than gmax in reality. 



Schwarzschild case, Ch can be evaluated using numeri- 
cal results of the greybody factor Q; e.g., Ch « 1-7 for 
d — 6. The greybody factors for spinning BHs are stud- 
ied for brane spacetimes [ol-fl^ whereas those for bulk 
gra vitons have never been completely evaluated (but see 
|13|). The general tendency is that the emission rate 
is enhanced significantly as the value of q or d is in- 
creased: The luminosity for q ^ 1 is ~ 100 times as 
large as that for q — (e.g., |11|). The angular momen- 
tum emission rate is more significantly enhanced as the 
spin is increased, indicating that a rapidly spinning BH 
spins down in a short timescale. Thus, strictly speaking, 
we should take into account two timescales: One is the 
evaporation timescale th and the other is the spin-down 
timescale referred to as ths- Reference [ll| shows that 
for q ~ 1, ths <C Th: ths for g = 1 is by a factor of ~ 100 
shorter than that for g ^ 1 , but th is not as short as ths ■ 
The ratio of th to tgw is 



TH 

Tgw 



Ch 



GW 



Mo 
mp 



(d-2)/(d-3) 



(74) 



and proportional to M^'^ 3) ^]-^gj.g ^ < [d—2)/{d — 

3) < 4/3 for d > 6. Thus, the ratio depends moder- 
ately on Mq for mp < Mo < lOmp. On the other hand, 
Ch/Cgw are in the range ~ 10"'^ - 1, depending strongly 
on the spin. This implies that the evolution of a mini BH 
formed after a particle collision depends strongly on the 
spin q (or equivalently, the impact parameter b). 

Now, we discuss the phenomenology of a mini BH as- 
suming that mp = 1 TeV and Mo = 10 TeV = lOmp, 
which are plausible values in TeV-gravity scenarios and in 
LHC. Note that the scenario is qualitatively unchanged 
as far as Mq < lO^mp. As a specific example, we further 
fix d = 6 in the following. For Mo = lOmp, tgw ~ 8tp 
and Th ~ 80tp tgw for q — 0- As the value of q 
is increased from q = 0, tgw becomes longer and th 
shorter (for a small value of g <C 1, ths is comparable to 
Th). Thus, two timescales, tgw and th, become identi- 
cal, Tgw = Th, at a value of g = qcq. The result in this 
paper indicates ^cq < gcrit, because tqw becomes very 
long (Cgw is of order lOO) for q ^ gcrit- 

Then, the phenomenology is classified into three types 



depending on the spin g; (i) O < 9 ^ ?eq, 



(ii) 9oq 



< 



q < 



gfrag, and (iii) q > qhag- Here, we assume a hypothetical 
critical value for the onset of fragmentation as (/frag > 1-6 
following the estimate in Sec. VI B, although it is not 
clear whether the fragmentation really occurs. Because 
of the high value of q, the fragmentation will proceed in 
a timescale of order Mq^^'^~^'^ ~ th ^ tgw, if it occurs. 

For the case (i), the standard hypothetical picture 
of BH evaporation holds; a formed BH emits gravita- 
tional waves and settles to a stationary state in a short 
timescale. Then, it will be evaporated by Hawking radi- 
ation. The formed BH is spinning, but not very rapidly, 
and the effect of the energy loss and spin-down by gravi- 
tational radiation reaction plays a minor role. The spin- 
down will proceed primarily by Hawking radiation, and 
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finally, a nearly nonspinning BH will be evaporated. This 
type of BH formation and evaporation will not be a dom- 
inant process in particle accelerators, because a mini BH 
with such a small spin {q < qcq < Qcnt) or with a small 
impact parameter will not be formed frequently; the for- 
mation rate will be by a factor of (gcq/'i'max)^ smaller 
than the total. 

For the case (ii), the BH will not be relaxed to a stable 
state by gravitational-wave emission, because of its long 
emission timescale. Such a nonstationary BH will start 
radiating quantum particles while emitting gravitational 
waves, and the signal of the quantum radiation in this 
phase is likely to be different from the idealized Hawk- 
ing radiation (which is the result only from a station- 
ary, axisymmetric BH). Because a BH with a large initial 
spin reaches a highly deformed state, non-Hawking-type 
quantum radiation will be enhanced for the higher spin. 
Obviously, an improved analysis for predicting the spec- 
trum of the quantum radiation is required. Assuming 
that the spin-down timescale by the quantum radiation 
is shorter than the evaporation timescale as in the Hawk- 
ing radiation 11[ , the BH spin will subsequently become 
smaller than qcrit after the substantial quantum radia- 
tion, and thus, the BH relaxes eventually to a stable, 
quasistationary, axisymmetric state and will stop emit- 
ting gravitational waves. After the condition q < qcrit is 
achieved, the quantum radiation process will be identical 
to the ordinary Hawking radiation. 

For the case (iii) , the timescale for the growth of bar- 
mode deformation and subsequent fragmentation is likely 



to be of order ~ M. 



l/(d-3) 




This may be as short as ths 



If fragmentation occurs, a rapidly spinning BH is likely to 
change to two slowly spinning BHs before the quantum 
(Hawking) radiation becomes the dominant dissipation 
process. A supportive evidence for this possibility is that 
the apparent horizon formed in particle collisions takes 
a peanut shape for b « 6max as shown in Figs. 5 and 6 
of [1^, which indicates that gravity combining the two 
particles is weak. If the fragmentation occurs (the BH 
horizon pinches off) because of quantum-gravity effect, 
two boosted BHs will be the outcome (binary will not be 
the result for d > 5). Thus, the Hawking radiation from 
the two boosted BHs may be observed subsequently, as 
two jets in particle accelerators. Here, we note that the 



two timescales, M, 



l/(d-3) 




and ths , may be as short as the 



Planck time rp for Mq w lOrnp, and thus, such a (semi) 
classical phenomena may be veiled by quantum-gravity 
effects in LHC. 



D. Issues for the future 

The present paper reports the results for the stability 
of rapidly spinning BHs against bar-mode deformation. 



We analyzed the BHs only with the spin q < 1.15 and 
did not study for the ultra spinning case with q ^ 1-5 
because it is not technically easy to perform a longterm 
simulation for such BHs. As discussed above, the evo- 
lution of the ultra spinning BHs may be qualitatively 
different from that for q I. Furthermore, such ultra 
spinning BHs may be the frequent outcomes, if the TeV- 
gravity hypothesis is correct. Clarifying the evolution of 
the ultra spinning BHs is obviously an important issue 
left for the future. 

The stability of black rings against bar-mode deforma- 
tion is one of the interesting issues. The analytic solu- 
tions for the black rings are found in five dimensions [13] , 
and numerical-relativity simulation in a similar manner 
to that in this paper may be possible. The black rings 
always have a high spin with q > gcrit for d = 5. Perhaps, 
they are also unstable against bar-mode deformation and 
evolve as a result of gravitational radiation reaction as 
far as q is not extremely large; if the spin is very large, 
fragmentation may occur. 

In the present paper, we start simulations in a highly 
idealized situation; we prepare nearly stationary, axisym- 
metric BHs and evolve them approximately in a quasis- 
tationary manner. In particle accelerators, however, the 
situation will be highly different. The rapidly spinning 
BHs after the particle collision will be highly nonstation- 
ary and nonaxisymmetric. Such BHs may evolve qualita- 
tively in a similar manner to that in the present analysis, 
but quantitative properties on the evolution process will 
be significantly different. To clarify the formation and 
evolution of such BHs, it is necessary to perform a sim- 
ulation started from high- velocity two-BH collision. Our 
goal is to successfully perform this simulation and to clar- 
ify the possible outcome in this setting. 
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